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Abstract 

Many disordered systems experience a transition from a fluid-like state to a solid-like state follow- 
ing a sudden arrest in dynamics called jamming. In contrast to jamming in spatially homogeneous 
systems, jamming in hoppers occurs under extremely inhomogeneous conditions as the gravity- 
driven flow of grains enclosed by rigid walls converges towards a small opening. In this work, we 
study velocity fluctuations in a collisional flow near jamming using event-driven simulations. The 
average flow in a hopper geometry is known to have strong gradients, especially near the walls and 
the orifice. We find, in addition, a spatially heterogeneous distribution of fluctuations, most strik- 
ing in the velocity autocorrelation relaxation times. At high flow rates, the flow at the center has 
lower kinetic temperatures and longer autocorrelation times than at the boundary. Remarkably, 
however, this trend reverses itself as the flow rate slows, with fluctuations relaxing more slowly 
at the boundaries though the kinetic temperatures remain high in that region. The slowing down 
of the dynamics is accompanied by increasing non-Gaussianity in the velocity distributions, which 
also have large spatial variations. 
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I. INTRODUCTION 



Granular flows can clog unpredictably, often with catastrophic consequences as seen in 
grain silo failures [I]. This phenomenon is an example of jamming, where a disordered system 
undergoes a sudden arrest in dynamics, leading to a transition from a fluid-like to a solid-like 
state [2]. Unlike the jams that develop in spatially homogeneous systems, however, jams in 
hopper or silo flows occur under extremely inhomogeneous conditions as a gravity-driven 
flow enclosed by rigid walls converges towards an orifice where stable arches can form[3H7]. 
Hopper flows are self-organized, being neither pressure nor density controlled, with a steady 
state that develops from an interplay between gravitational forcing, constraints due to the 
walls, and the interactions between grainsjE]. The flow is also driven by body forces, and the 
boundary is thus not a source of energy. Instead, the boundary imposes constraints on the 
flow volume, is a source of frictional interaction, and randomizes the flow through grain-wall 
collisions. The dynamical principles that lead to jamming in such self-regulated flows are 
not well understood. 

Experimental investigations of two-dimensional (2D) hopper flows show that these flows 
remain collisional near jamming[9j, and show large force fluctuations [10J and transient arch 
formation |5j. In previous studies [TT1 [12], we have shown that collisional flows with no ex- 
tended contacts develop transient force chains that are are sustained by correlated collisions. 
The origin of these correlations is the inelasticity of the collisions, and the stress is purely 
due to momentum transfer. Focusing on the homogeneous flow in the bulk of the hopper, 
we established the emergence of growing length and time scales as the flow approached 
jamming IT^] . Flows in hoppers, however, are characterized by strong gradients, such as 
a shear-layer at the walls [T5|IT6]. and rapid flow regions near the orifice. In addition to these 
heterogeneities in the average flow profile, the fluctuations in the flow have well-defined spa- 
tial structure, such as transient arches jS]- In this work, we focus on the dynamics in regions 
with strong gradients and their vital role in controlling the approach to jamming. 

Using event-driven simulations of a purely collisional, gravity-driven flow, we monitor the 
temporal variation of the flow in boxes that are a few grain diameters in size. We observe 
that the approach to jamming in the collisional non-equilibrium steady state (NESS) is 
signaled by the flow becoming intermittent with time intervals spanning many collisions in 
which the flow alternately slows down and speeds up. During some of these intervals, the 
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flow is slow (fast) throughout the hopper but in other instances, the slowing down (speeding 
up) is localized in the bulk or near the boundaries. These fluctuations are particularly strong 
in hoppers with small openings f9]. In addition to the marked gradients in the flow field, the 
fluctuation-statistics of the velocity show large variations both along and transverse to the 
flow direction. More importantly, the spatial patterns evolve with flow rate and can change 
quite dramatically as the flow approaches jamming. In particular, we observe a flowing 
regime, and a pre-jamming regime distinguished by the spatial behavior of the velocity 
relaxation. In a rare jamming event that spans millions of collisions, we see evidence of the 
boundaries frustrating the tendency to develop a uniform, steady flow in response to the 
gravitational forcing. 

II. DESCRIPTION OF SIMULATION 

Our event-driven simulation is based on a particle dynamics model used by Denniston 
and Li [T7j and is described in detail in earlier papers [TT1 IT2] . The system consists of 1,000 
non-deformable disks falling under gravity in a 2D rectangular hopper with a tapered base 
(see Fig. [I]). The disks are bidisperse, with diameters d and 1. 2d; collisions between particles 
are inelastic but frictionless, so momentum transfer occurs along the center-to-center vector 
of colliding particles. The transformation in relative velocity between colliding particles i 
and j is defined in terms of a coefficient of restitution /a: 



where u , u' are the particle velocities before and after the collision, and q is a unit vector 
along the center-to-center direction. All inter-particle collisions become elastic below a 
relative velocity threshold u cut to avoid inelastic collapse [18]. Particle- wall collisions are 
inelastic, with frictional drag at the walls modeled by a tangential coefficient of restitution 
// W aii, allowing the flow of grains through the outlet at the bottom of the system to reach 
a steady-state. Particles exiting the system at the base are reintroduced at random lateral 
positions at the top of the hopper. 

All collisions in this event-driven simulation are instantaneous, and particles cannot form 
temporally extended contacts with one another. The simulation moves forward from collision 
to collision, and instantaneous snapshots of the system are reconstructed from the state of the 
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system before and after the relevant collisions to determine particle positions and velocities 
at equally spaced time intervals. The overall flow rate is set by the particle-wall friction 
coefficient and the size of the opening and varies linearly with opening size, as we show in 
Figure [T} The mass of the smaller grains and the acceleration of gravity g are set to 1, 
and all lengths are expressed in units of the smaller particle diameter d. In these units, the 
length and width of the rectangular region are 76.5 and 20 respectively, the angle of the 
tapered base is 45 degrees, and the size of the opening at the base is varied between 10 and 
4.5. The simulation parameters for the collisions are ji = 0.8, /i wa n = 0.5, and u cut = 10~ 3 . 
The time taken for free fall of a particle through its own diameter is a/2, and the average 
time between collisions for a given particle is on the order of 10~ 3 . 

The development of extended contacts is, of course, important for sustaining a jam. It 
is not clear, however, if the approach to jamming is qualitatively altered by the presence 
of extended contacts. Frequent collisions can mimic extended contacts[TlJ [12], and one 
can construct an effective potential for hard spheres by coarse-graining configurations over 
time intervals large compared to times between collisions[19, 20j. In this paper, we discuss 
signatures of an approaching jam in collisional flows, assuming that a suitable coarse-graining 
can map these on to flows that develop extended contacts. The majority of our results are 
from the intermittent flow regime where we observe clearly defined intervals of slow and fast 
flows. In addition, we have studied one jam that lasts for 10 7 collisions, and characterized 
the structure of the flow during its creation and breakup. A permanent jam cannot be 
created in our model of collisional flow since we prevent inelastic collapse [TS] . 



III. RESULTS 



To study fluctuations, we construct a velocity field, v(r,t) by coarse graining individual 
particle velocities ui over a box of size 2d x 2d: 

1 N 

Hr,t) = -Y,^(t) (2) 

i=i 

where the sum is over the N particles whose centers lie inside a box centered at r = (x, y) at 
time t. The entire hopper is subdivided into non-overlapping boxes, starting at the base of 
the hopper and proceeding vertically upwards. We employ two different types of averaging to 
measure fluctuations. A time-averaged flow-field, V{r) = J dtv(r, t) may be used to measure 
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FIG. 1: (Color Online) Left: Spatial contour map of V y (r), illustrating a gradient-free region in 
the center, a shear layer bordering the walls, and a region of rapid flow near the orifice. Center: 
A snapshot of all the particle positions at a given instant in time, showing fast particles with 
Ui,y(t) > V y (t) (red) and slow particles with Ui >y < V y (t) (blue). Upper right: Mean flow rate, 
% w = J dtV y (t), measured in simulation units, shown as a function of the outlet size measured 
in units of d. Lower right: Locations ({x,y}) of three boxes that exemplify three very different 
regions of the hopper. 

temporal fluctuations in each box; and a spatially averaged flow-field, V(t) = J drv(r, t) may 
be used to measure instantaneous variations of the velocity. Both measures are illustrated 
in Figure [TJ which shows the spatio-temporal heterogeneity of the velocity field. In most of 
the analysis we describe in this paper, we use the box-specific definition, V(r). 

The leftmost panel of Figure [I] shows the spatial pattern in V y (r). Particles enter the 
hopper at the top and large gradients are seen in the direction of the flow. This region 
transitions into the "bulk" where V y (r) is gradient-free in the direction of the flow, but a 
marked interface with a large gradient separates the bulk flow from the shear layer along 
the vertical wall. This bulk flow gives way to a flow with strong gradients, transverse to 
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and along the flow direction, close to the orifice. The length scales associated with the 
gradients of V y (r) are not well understood. The width of the shear layer has been related 
to a stress-activated mechanism [T5], and we are currently investigating the relation between 
stress fluctuations and velocity gradients in our simulations [2l]. 

The center panel of Figure [l] shows an example of instantaneous variations in the velocity 
field. Tracking the velocities of individual particles, Ui(t), we label them as fast (ui tV (t) > 
V y (t)) and slow (u ity (t) < V y {t)). The snapshot shows large, vertically-extended clusters of 
fast moving particles, demonstrating non-trivial spatio-temporal correlations in the flow. 

The average flow field, V y {r) identifies three interesting regions of the flow, if we set aside 
the region at the top of the hopper where particles are entering. To analyze fluctuations 
in these regions, we choose three representative boxes: a box at the orifice, a box in the 
bulk flow region, and one in the shear layer; these three are shown in the lower right panel 
of Figure [TJ In addition, to better characterize the spatial variation of the fluctuations, 
we study their statistics in boxes along the central vertical column, and a horizontal row 
running through the bulk flow region. These two cuts span the regions of strong gradients in 
V y (r). The mean flow rate, vr ow = J dtV y (t), is controlled by the size of the hopper opening, 
and as the upper right panel of Figure [T] shows, varies linearly with the outlet width. 

A. Velocity Autocorrelations 

The approach to jamming is marked by large temporal variations in v y (r, t). As we show 
in this section, the timescales associated with velocity fluctuations provide the strongest 
signal of an incipient jam. We study the statistics of the temporal fluctuations through the 
velocity autocorrelation function, C(r,t) = (v y (r,t)v y (r,0)). As an example of the behavior 
of C(r,t), the upper panel of Figure [2] shows its variation with flow rate in the shear- layer 
as measured in the box closest to the vertical wall. If we measure the autocorrelation time 
by C(r,r(r)) = 0.1, then r in this box increases by a factor of 2 as the outlet size changes 
from 10 to 4.5. The spatial variation of the relaxation times of the velocity fluctuations 
is represented by the field T(r), which is shown in Figure [2] for the fastest and the slowest 
flow rates. The plots are symmetrized about the central vertical line, so there is meaningful 
information only in half of each plot. 

The lower left panel of Figure [2] illustrates the pattern of r(r) at an opening of 10.0. 
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This fast-flowing NESS is characterized by regions of fast relaxation close to the boundary, 
and increasingly slow relaxations as one moves towards the center of the hopper. At slower 
flow rates, a clear reversal in trend is evident from the lower right panel of Figure [2] At 
opening size 4.5, r(r) is longest at the walls, and shortest at the center. The change in 
the pattern of relaxation times signals a qualitative change in the behavior of the flow as it 
approaches jamming, and is the most significant dynamical signature that we have observed 
in our simulations. 

The near-jamming pattern poses a puzzle if we view the NESS as a fluid characterized by 
a granular temperature [22J. The kinetic granular temperature is proportional to the variance 
of the velocity distribution. As shown in our earlier paper [13], and detailed later in Table |ll[ 
the kinetic temperature increases from the center to the shear layer with a sharp peak 
very close to the vertical wall. This trend does not change with flow rate. Near jamming, 
therefore, velocity fluctuations relax more slowly in regions with higher kinetic temperature. 
This behavior is at odds with our intuition, based on thermal fluids, of slower relaxations at 
lower temperatures. For the fast flows, the hopper-NESS does show the expected behavior. 
We conclude that there is a pre-jamming regime in collisional hopper flows where the NESS 
behaves neither like a fluid nor a solid. 

If we assume that the Green-Kubo relations |23j hold for the NESS in the hopper, then 
the autocorrelation time is related to the diffusion constant 1/r = D oc ksT/r), where T is 
the kinetic temperature and r\ the viscosity. The viscosity is related to fluctuations in the 
stress through another Green-Kubo relation. At the fastest flow rate, we observe that D is 
small at the center and increases towards the boundaries, consistent with the trend in T. 
The fast-flowing NESS thus qualitatively resembles a fluid with temperature gradients. For 
slow flows, D is smaller at the boundaries than the bulk. One possible explanation is an 
anomalously large viscosity rj near the boundaries, which just trades one puzzle for another 
since viscosity should also decrease with increasing temperature. We do not understand 
the origin of the slow-relaxation pattern close to jamming. A reasonable hypothesis is 
that boundaries frustrate the downward motion of the gravity-driven flow leading to slow 
relaxations. In future work, we will study stress fluctuations [21] and their relationship to 
the velocity fluctuations in order to better understand the pre-jamming NESS. 

The pattern of relaxation times indicates that jamming is induced by the boundaries and 
propagates into the bulk. As the flow slows, we also observe an increasing frequency of events 
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FIG. 2: (Color Online) Upper panel: C(r, t) measured at the boundary box at different outlet 
openings. Lower panel: Spatial variation of r(r) at outlet size 10.0 (left), which was the fastest 
flow rate simulated, and at outlet size 4.5 (right), which was the slowest flow rate simulated. 

with zero particles at the outlet. These events are signaled by a peak at zero velocity in the 
velocity distribution function, seen in the upper right panel of Figure [3j The percentage of 
these "zero" events remains small, varying from 0.14% at opening 10.0 to 0.67% at 4.5, as 
might be inferred from the time-trace of the velocity on the upper left side of this figure. 
The boxes in the shear layer also show a large number of zero and negative velocity events, 
as shown in the two lower panels of the same figure, consistent with the picture of a jam 
originating at the boundary. The shape of the distributions in these two boxes are distinct 
and will be discussed in Section C. In a future publication, we will analyze the negative 
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FIG. 3: Illustrative examples of the time series of the velocity field, v y (r,t) at opening size 4.5. 
(a) Time series measured in the box at the outlet, (b) the corresponding distribution; (c) Time 
series measured in the box at the vertical wall, (d) the corresponding distribution. 



fluctuations from the perspective of non-equilibrium fluctuation relations 



B. First Passage Time Distributions 

In the pre-jamming regime, our analysis of the autocorrelation functions paints a picture 
of a collisional flow with fluctuations that do not fall within any familiar rubric of fluid 
flow. In this section, we analyze first passage processes in which the velocity transitions 
from below to above the average velocity, and those that go from above to below, to gain 
more insight into the nature of the intermittent flow that develops very close to jamming. 
It is well known [25J that for random processes, autocorrelation functions and first-passage 
probabilities provide complementary information. 

To study the first passage probabilities, we introduce a clipped variable at each box, 
defined as s(t) = v/\v\, which changes sign at the zeroes of a scaled velocity variable v. This 
scaled velocity is defined for each box as v = (v y — (v y ))/a where (v y ) = V y {r = XboxiVbox) 
and a is the standard deviation obtained from the time series, v y ((r = Xbox,ybox),t). We 
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TABLE I: This table shows, the average first passage times for the velocity at a particular location 
to go from above to below the mean ((r+)), or or from below to above the mean ((r_)). The rows 
correspond to different flow rates, and each pair of columns corresponds to the average times at 
the outlet, boundary, and bulk boxes. 

then define P+(t) (P_(£)) as the probability that this function does not change sign during a 
time interval of length t, i.e. retains a value s(t) = +1 (s(t) = —1) and contains no zeroes. In 
other words, P+(t) (P_(i)) measure the probability that a fluctuation that is faster (slower) 
than the mean, lasts for precisely a time interval t. We examine these processes in the three 
representative boxes, highlighted in Figure [TJ which represent the different regimes of V y (r), 
and where we also look at the distributions of the velocity and its autocorrelation function. 
These distributions, shown in Figure |4j appear as power laws cutoff by exponentials for all 
flow rates, and in all the boxes studied. The cutoff times increase as the flow rate decreases, 
implying that the average first passage time increases as the flow slows down, whether we 
are looking at the time characterizing the switch from above to below the mean ((r + )), or or 
from below to above the mean ((r_)). The times (r + ) and (r_) obtained from exponential 
fits to the tails of the data are shown in Table |T} The distributions unambiguously show an 
increasingly intermittent flow. We speculate that the trend in average first passage times 
indicates some type of localization in trajectory space, akin to proposed mechanisms for the 
glass transition |26j. 
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FIG. 4: (Color Online) Distributions of first passage times, P+(t) and P—(t): Outlet Box (top 
panel: a, b) and Boundary box (bottom panel: c, d). The cutoff times obtained from fits to the 
tails of these distributions are shown in Table HI 

C. Velocity PDFs 

The probability distributions functions (PDFs) of the vertical components of the velocity 
become increasingly non-Gaussian as the flow reaches the pre-jamming regime. A regime of 
non-Gaussian fluctuations has been observed in experiments on two-dimensional silos below 
a certain orifice size[27]. In this section, we analyze the spatial variation of these PDFs. 
We have constructed box-specific distributions of v y (r,t) and of the scaled velocity variable 
v(r,t) = (v y (r,t) — (v y ))/a by monitoring the time evolution of the velocity in each box. 
Figure [5] shows distributions of the scaled and unsealed velocities for a column of boxes 
parallel to the flow and including the orifice, and a row of boxes transverse to the flow 
and including the shear layer. All of the distributions exhibit significant non-Gaussianity, 
measured quantitatively by the skewness and kurtosis. Table [IT] summarizes the values of 
the mean, the standard deviation and skewness for these boxes. 

Different classes of velocity distributions in different regions have been noted in event- 
driven simulations of flow in three-dimensional silos that are not tapered [25]. We find marked 
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FIG. 5: (Color Online) Velocity Distributions at the smallest outlet size, 4.5d. Top Panel: (Left) 
The distributions of the scaled vertical velocity v = (v y — {v y ))/a for the column of boxes shown 
in the inset; (Right) distributions of the unsealed velocity v y for the same set of boxes. Bottom 
Panel: (Left) distributions of v for the row of boxes shown in the inset; (Right) distributions of the 
unsealed v y for the same set of boxes. The corresponding values of the mean, standard deviation 
and skewness are summarized in Table UH 

spatial variation of the velocity fluctuations near the orifice and the vertical walls, whereas 
the three-dimensional flows showed significant variations only near the top of the hopper. It 
is not clear whether these differences stem from dimensionality, the geometry of the hopper 
or both. Our results are relevant for experiments performed in a 2D tapered geometry [9]. 

Figures [5] and [6j which show velocity distributions at different flow rates, illustrate some 
remarkable features of the velocity distribution. At all flow rates, and in all regions that 
we have examined, except the shear layer spanning 3 — 4 grain diameters bordering the 
vertical wall, the distributions of the vertical component of the velocity are described re- 
markably well by a universal form, the generalized Gumbel (GG) distribution [29]. The GG 
is a generalization of distributions found in extreme value statistics [29] that appears to be 
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FIG. 6: (Color Online) Comparison of the velocity distributions at different flow rates, Vfi ow . Top: 
Distributions of the scaled velocity, v (left), and the unsealed velocity, v y measured at the outlet 
box. Bottom: Distributions of the scaled velocity, v (left), and the unsealed velocity, v y measured 
at the vertical boundary box. 

more broadly applicable: it has been observed in a simulated low density granular gas of 
hard disks j3(]l El] where it describes the energy distribution, and in sheared systems close to 
jamming where it describes entropy fluctuations |32j. The GG distribution is entirely deter- 
mined by moments of the measured distribution, namely the mean, skewness and standard 
deviation: 

IL(y) = Ke a[x - eX \ with x = -b{v + z) . (3) 

Here K, b, z are all functions of the parameter a, which is in turn related to the skewness 
of the distribution: \(v 3 )\ = 1/ ' \fa. The functions appearing in Equation [3] depend only on 
a: 

Id 2 lnT(a) 1 / , . dlnT(a)\ ba a 

where T(a) is the Gamma function. 

It should be emphasized that the GG is not a fitting form, but is simulated from measured 
properties, and compared to the measured distribution. As an example, the figure in the 
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TABLE II: The mean, (v y ), standard deviation, a = \J (v 2 ), and skewness, (£ 3 ) of the velocity 
distributions shown in Figure [5] The panel on the right compares the distribution of v at the 
orifice to the best Gaussian fit (blue, dashed) and the simulated GG distribution (red, solid). 

right panel of Table [H| shows a comparison between the best Gaussian fit (blue, dashed), 
the GG (red, solid), and the measured distributions of v at the outlet. At a given flow rate, 
the parameter a characterizing the GG, is smallest closest to the boundaries and increases 
towards the center of the hopper. The distributions are, therefore closer to a Gaussian at 
the center where the flow rate is highest. The skewness also changes sign as a function of 
vertical position of the box, being negative and large at the orifice and positive but smaller 
in the central region. A negative skewness indicates a larger weight in the distribution below 
the peak than above it. At a given position in the hopper, the GG parameter decreases with 
decreasing outlet width. The GG form describing the measured velocity distributions are, 
therefore, farthest from a Gaussian at the slowest flow rate and closest to the orifice. The 
GG parameter a provides a quantitative measure of the variation of the fluctuation-statistics 
with flow rate. 

The appearance of a universal, non-Gaussian distribution suggests that there is an un- 
derlying physical principle leading to this particular form of non-Gaussian behavior. We do 
not yet have any insight into this physical mechanism, however, we expect that comparing 
to other near-jammed systems will be a productive avenue to pursue. 

The shape of the distribution within the shear layer at the vertical wall, while strikingly 
non-Gaussian, cannot be described by the GG form either. There are significant negative 
velocity events that signal particles moving opposite to the direction of the driving, and the 
distribution of these negative events is very close to exponential. The edge of this exponential 
behavior seems to be pinned at zero velocity for a range of outlet openings, as seen in the 
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lower panel of Figure |6j The appearance of significant negative-velocity events seems to 
be a hallmark of near-jammed systems [321 Intriguingly, there is remarkable similarity 
between the velocity distribution in the shear layer of our model hopper flow and that of 
a polar granular rod moving in a vibrated bed of granular spheres [33]. Negative events 
are also observed at the orifice of the hopper, however, their contribution to the velocity 
distribution appears to be consistent with the shape of the GG, as shown in the right panel 
of Table [TTJ Zero flow events at the outlet box, not shown in Table (TTJ, appear as a narrow 
peak superposed on top of the GG distribution as seen earlier in Figure [3] these occur when 
there are no particles in the outlet box in a given time interval. 

In the next section, we describe our analysis of a jam that persists for 10 7 collisions, 
and provide evidence for the appearance of spatial structures with growing length scales, 
including vortices. 

D. Jam Event 

So far, our discussion has focused on fast flows and the pre-jamming flow which is an 
interesting NESS that cannot be cleanly labeled as a fluid or a solid. In this section, we 
analyze a jam where the the flow stops for a time interval spanning 10 7 collisions. This event 
offered us the rare opportunity of studying a NESS that is as close to a solid as possible in 
purely collisional flows. 

We observed that the jam was due to the formation of a classic arch at the outlet. We 
tracked the particles that ended up in the arch, Figure [7], and found that arches are formed 
at the outlet and not pre-formed upstream in the flow. In the simulation, as in experiments 
on photoelastic beads[3], arch-forming particles originate in disparate regions of the system. 
These similarities with experiments suggest that a jam in a collisional flow might not be 
that different from one in a physical system with extended contacts, and encouraged us to 
analyze the anatomy of this jam. We studied two aspects of the jam: (1) appearance of 
vortices in the velocity field and (2) geometry of spatio-temporal clusters. 

We tracked the time evolution of the velocity field v through the formation and disso- 
lution of the jam. The upper panel of Figure [8] shows the velocity field at two instants in 
time leading up to the jam, and clearly demonstrates the increasing vorticity of the flow. 
Comparing the snapshot on the left to the one on the right, which is at a later time, we 
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FIG. 7: (Left) (Color Online) A snapshot taken approximately 10 3 collisions before the arch 
formation, showing the position of the particles that end up forming the arch (Right) responsible 
for the jam. Tracking the motion of these particles, there motion is found to bear a striking 
resemblance to the process of arch formation seen in experiments |5j. 

see that vortices nucleate at the edges of the hopper, and propagate into the bulk. We are 
currently in the process of measuring this behavior systematically [33]. To our knowledge, 
vorticity of self-organized flows such as the hopper flows have not been investigated much 
nor have they been correlated with jamming. In contrast, length scales associated with 
vorticity have been studied in chute flows [35], and have been related to jamming and flow in 
that geometry through a scaling analysis. Our results suggest that vorticity could be related 
to the dynamical arrest of hopper flows. The negative velocity events that we observed in 
the intermittent flow regime can also be related to vortices, but are transient, and unlike 
the ones in Figure [8j do not have any significant spatial extent. 

We have also analyzed the spatial extent of 'slow clusters'. We define slow boxes as ones 
that have, at a given instant, a vertical velocity smaller than half of the average flow rate 
in the box. Clusters are identified by nearest-neighbor connectivity and then tracked as a 
function of time. The lower panel of Figure [8] shows two snapshots of cluster formation in 
the time leading up to the jam. Each of these cluster images corresponds to the same instant 
in time as the vorticity image in the panel immediately above. The two images on the left 
demonstrate that there is vorticity associated with the regions where we see large clusters. 
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FIG. 8: (Color Online) The upper panel of this figure shows the instantaneous velocity field in 
the lower region of the hopper at two times leading up to the jam. The lower panel corresponds 
to the same two instants in time, and highlights boxes in which the instantaneous vertical velocity 
is less than half the mean in that box. The figures on the right are closer in time to the jam. 

The images on the right are closer to the instant of the jam, and the cluster spans nearly the 
entire region. To obtain quantitative information, we plot cluster function of time 

both during the formation and the breakup of the jam. Figure [9] shows that formation and 
breakup of the jam is signaled by a rapid change in cluster area. This is not surprising and 
corroborates the picture that a stable arch at the outlet supports the weight of the grains 
above while large regions of slow flow develop. 

The clusters that are observed during the jamming event are an example of a dynamically 
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Time before jam Time after jam 

FIG. 9: Scatter plot of the area covered by connected clusters, defined by boxes in which the 
flow is less than half the average, plotted against time before (on the left) and after (on the right) 
the jamming event. The duration of the jam that separates the two figures is about 24 simulation 
times. 

generated structure. Similar structures have been seen in experiments [9], so we do not think 
this is an artifact of our simulation. In our collisional model, these slow flow regions are 
nowhere near static, instead, the grains continually undergo collisions. The clusters retain 
their identity while their constituents participate in a highly correlated yet rapid dynamics. 
Our ultimate aim is to understand the dynamical principles that lead to the emergence of 
these structures, and their stability. In the next phase of this work, we will investigate the 
formation and breakup of stable, long-lasting jams by generating many realizations of our 
simulation at the smallest opening[34j, and by using molecular dynamics simulations with 
particles interacting via Hertzian contacts [36J. 
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IV. CONCLUSION 



We have presented an analysis of the jamming transition in purely collisional 2D hopper 
flows. In our simulations, we identify new dynamical signatures of the approach to jamming 
in this self-regulated flow. In particular, we analyze features that are strongly influenced 
by the boundaries, and find a striking change in the nature of the flow as the opening size 
is reduced. The spatial distribution of velocity-autocorrelation times shows that, in a pre- 
jamming regime, the flow has a distinctly heterogeneous character. There is a glassy layer at 
the walls, characterized by anomalously slow velocity relaxations, coexisting with a central 
flowing region. This is in contrast to higher flow rates, where the flow resembles a thermal 
fluid with strong spatial gradients. The flow is strongly intermittent in the pre-jamming 
regime, with well-defined periods of slow and fast flow. From measurements of first-passage 
time distributions, we deduce that velocity fluctuations that dip below the mean persist 
longer as jamming is approached. The effect is strongest at the shear layer, where we also 
see the development of a glassy region. 

We also observe marked vorticity in the flow as an extended jam forms and then breaks up. 
The vortices nucleate at the corners of the hopper and extend inwards, ultimately causing 
complete arrest. The vorticity is spatially correlated with the solid-like regions, suggesting 
that the slow relaxation of the velocity is related to the vorticity. In our earlier work|12j. we 
saw clear evidence of stress chains in this purely collisional flow. These stress chains were 
predominantly anchored on the side walls, and it is plausible that they are being supported 
by the glassy layer. In future work we will investigate the mechanism by which the more 
glassy boundary layer sheds vortices into the free- flowing bulk layer to generate jams. One 
question that arises is whether extended contacts are necessary for the onset of jamming; 
our results here suggest otherwise. We will carry out molecular dynamics simulations and 
compare them to the event-driven results in order to fully answer this question. 
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